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Abstract 

The adiabatic inhomogeneities of the scalar curvature lead to a compressible flow af- 
fecting the dynamics of the hydromagnetic nonlinearities. The influence of the plasma on 
the evolution of a putative magnetic field is explored with the aim of obtaining an effective 
description valid for sufficiently large scales. The bulk velocity of the plasma, computed in 
the framework of the ACDM scenario, feeds back into the evolution of the magnetic power 
spectra leading to a (nonlocal) master equation valid in Fourier space and similar to the 
ones discussed in the context of wave turbulence. Conversely, in physical space, the mag- 
netic power spectra obey a Schrodinger-like equation whose effective potential depends on 
the large-scale curvature perturbations. Explicit solutions are presented both in physical 
space and in Fourier space. It is argued that curvature inhomogeneities, compatible with 
the WMAP 7yr data, shift to lower wavenumbers the magnetic diffusivity scale. 
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1 Introduction 



The evolution equations of the magnetic power spectra in a conducting fluid have been the 
subject of extended studies with particular attention to the regime of high Reynolds numbers 
(see, e.g. [H El El II])- In standard treatments the flow is often passive (in the sense that 
it is given as an external input of the problem), it is isotropic and and it is assumed to be 
incompressible in various situations ranging from the usual one-fluid description provided 
by magnetohydrodynamics (see, e.g. [4J) to the apphcations of gyrotropic turbulence [3]. 
Exceptions to the previous statement are, for instance, the compressible turbulence [5] or 
the hydromagnetic evolution in the presence of acoustic disturbances. 

Prior to photon decoupling, the physical properties of the primeval plasma can be directly 
scrutinized by means of the measurements of the Cosmic Microwave Background (CMB) 
temperature and polarization anisotropics. These observations show, with an accuracy which 
is sufficient for the present considerations, the predominance of the adiabatic curvature 
perturbations over any other (possibly subleading) entropic contribution [HI E] (see also 
[8] for a recent review). Indeed, the results of the WMAP observations (see, e.g. the 
WMAP 7yr data release [9l [TOl [11]) imply that the position of the first acoustic peak of 
the temperature autocorrelations and the position of the first anticorrelation peak in the 
temperature-polarization power spectra are in a fixed ratio (i.e. approximately 3/4). Such 
a ratio can be easily derived by assuming a dominant adiabatic component of curvature 
inhomogeneities in the initial conditions of the Einstein-Boltzmann hierarchy. Tensor and 
vector modes play a negligible role (see [121 US] and references therein) but are anyway absent 
in the vanilla ACDM scenario (where A stands for the dark energy component and CDM 
for the cold dark matter component). The adiabatic nature of the large-scale curvature 
perturbations determines also the statistical properties of the fiuid fiow: the baryon-photon 
and the cold dark matter velocities are Gaussian, isotropic and irrotational. The bulk velocity 
of the plasma (i.e. Vh in what follows) coincides, at early times, with the baryon velocity [E] 
(see also [U]). When photons and baryons are tightly coupled the photon-baryon fiow is also 
compressible and the typical amplitude of the corresponding velocity field is determined by 
the large-scale curvature perturbations. Finally the correlation functions of the bulk velocity 
of the plasma is exponentially suppressed for large times beyond a typical time scale, denoted 
by Td and related to the diffusive (or Silk) damping. 

The main motivation of this paper can be be summarized as follows. Prior to decoupling 
the Prandtl number (i.e. Pr^s,gn = -Rmagn/-Rkin) is -P'^magn = C(10^°) with i?kin < 1 and 
-Rmagn ^1. In such a system the fluid flow is computable (for instance in one of the 
popular versions of the ACDM scenario) and, in some sense, even accessible experimentally. 
Consequently, a natural question to ask is wether it is possible to account for the effect 
of the large-scale flow on the evolution of the magnetic power spectra. This problem is 
formally analog to (but physically very different from) some of the themes mentioned in 
the flrst paragraph of this introductory section. In the recent past close attention has been 
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paid to the analysis of the effects of large-scale magnetic fields on the evolution of curvature 
perturbations and on the CMB anisotropies (see, e.g., [HI [16] and references therein). In 
this analysis the complementary (and to some extent inverse) point of view will be taken: 
given the properties of the large-scale flow in the framework of the ACDM paradigm we 
ought to know how is the evolution of a putative magnetic field modified, in particular for 
large-scales. 

A similar problem arises in the statistical treatment of conducting fluids [H l2l E] where an 
interesting approach has been pioneered by Kazantsev [17\, Kraichnan and Nagarajan [18] 
(see also [HI [20]) and Vainshtein [21]. The common aspect of the approaches of [171 [13 [21] 
is that the flow is assumed Gaussian, isotropic and incompressible: these are also the three 
main physical analogies with the situation addressed in this paper. 

In spite of these relevant analogies there also two important differences. The first one 
is that the fluctuations of the geometry (as well as the effects related to the expansion of 
the background) do not play any role in the hydromagnetic treatment of conducting fluids; 
conversely, prior to decoupling, the large-scale flow is exactly determined by the fluctuations 
of the spatial curvature. In this sense the flow discussed in the present analysis is not passive 
but rather computed from the large-scale curvature perturbations. The second physical 
difference concerns the hierarchy between the Reynolds numbers which are simultaneously 
high (i.e. i?kin ^ 1 and -Rmagn ^ 1) in the standard hydromagnetic treatment^ p/Tj (see 
also [TSl 121] )• This is not the case after electron-positron annihilation and the smallness of 
i?kin is partly related to the properties of the bulk flow whose evolution is described within 
the same perturbative treatment used for the large-scale curvature fluctuations. 

The layout of the present paper is the following. In section[2]the set of governing equations 
will be introduced in its general form. In section [2] the ACDM fluid flow will be fed back 
into the evolution of the magnetic fields leading to an explicit integrodifferential equation 
valid for a compressible flow. The evolution of the magnetic power spectra is derived both 
in Fourier space and in physical space (section [4]). Explicit solutions are derived in section [5] 
Section [6] contains the concluding remarks. In the appendixes useful technical results have 
been collected for the interested readers. 

2 Governing equations 

The hydromagnetic equations in the presence of the fluctuations of the geometry have been 
deduced in different gauges (see, e.g. [HI [16] and references therein). Here the essentials of 
the curved space description of the predecoupling plasma will be briefly reviewed. The ap- 
proach of this paper is, in some sense, inverse of the analysis of [HI [16]: we are here concerned 

^ In spite of the fact that the three approaches of [HI [21] are conceptually related, the analysis of 
Kazantsev [T7] seems more appropriate for the present ends and will therefore constitute the basis for the 
generalizations developed in the present paper. 



3 



with the problem of deducing an averaged description for the evolution of the magnetic power 
spectra when the flow is determined in the framework of the ACDM paradigm. 



2.1 Curvature perturbations 

In a conformally flat background metric = a?{T)ri^y the relativistic fluctuations of the 
geometry are given by 

5s9m{x, t) = 2 a^(r)0(f , r), SsQijix, r) = 2 a^(r)V^(f , T)Sij, (2.1) 



in the longitudinal coordinate system. In Eq. (12. ip 6s denotes a metric perturbation which 
preserves the scalar nature of the fluctuation since, in the ACDM paradigm, the dominant 
source of inhomogeneity comes from the scalar modes of the geometry while the tensor 
and the vector inhomogeneities are absent. The scalar inhomogeneities are customarily 
parametrized either in terms of the curvature flutuations on comoving orthogonal hypersur- 
faces (conventionally denoted with TZ) or in terms of the curvature perturbations on uniform 
density hypersurfaces (conventionally denoted by (). Both TZ and ( can be defined in terms 
of the variables and ip appearing in Eq. (12.11) : 

n{n<j> + d^^) v^^ .... 

where "H = aH and H is the usual Hubble rate. The variables TZ and ( are gauge-invariant 
and the second relation of Eq. (12. 2p stems from the Hamiltonian constraint written in the 
longitudinal gauge of Eq. (12. ip . For a more thorough definition of large-scale curvature 
perturbations in different coordinate systems see fi5\ [23| I2l] and references therein. The 
large-scale curvature perturbations are then assigned in terms of the two-point function of 
7^=K, in Fourier space and prior to matter-radiation equality where curvature perturbations 
are approximately constant, at least in the case of the ACDM paradigm: 

In Eq. (12. 3p An = 2.43 x 10^^ is the amplitude of the power spectrum at the pivot scale 
kp = 0.002 Mpc~^ and Ug = 0.963 is, by definition, the scalar spectral index. The numerical 
values of An and Ug are determined from the WMAP 7yr data analyzed in the light of the 
vanilla ACDM scenario; the values of the remaining parameters are given by [H [TOl ITT] : 

{Q^o, Qco,^deo, ho, n^, ere) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088), (2.4) 

where, as usual, Qxo denotes the present critical fraction of the species x; ho is the Hubble 
constant in units of 100 km/sec/Mpc and ere is the optical depth to reionization. 
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2.2 Two-fluid equations 

The comoving electromagnetic fields and the comoving concentrations of electrons and ions 
are given by 

E{x, t) = a\T)S{x, r), B{x, r) = a'^{T)B{x, r), 

ni{x,T) = a^{T)hi{x,T), ne{x,T) = a'^{T)he{x,T). (2.5) 

Thus Maxwell's equations read 

V ■ E = A7ce{ni-n^), V ■ 5 = 0, (2.6) 

V X E = -drB, V X B = 47re(ni Vi - v^) + drE. (2.7) 

The velocities of the electrons, ions and photons obey, respectively, the following set of 
coupled equations: 

d^v, + nve = "-j[E + {Te X fi] - V0 + -^aT^,{v^ - v,) + aT,i{vi - v,), (2.8) 

Pea* 3 Pe 

drVi + nvi = — j[E + v,xB]-V(l) + -^ar^i(w, - v;) + aTei-ive - v^, (2.9) 
Pi a o Pi Pi 



drV^i = — -V5^ — V0 + ar^i('yi — Uy) + aT^eive — v^)- (2-10) 

In Eqs. fl2.8p -( l2TT0|) the relativistic fluctuations of the geometry are included from the very 
beginning in terms of the longitudinal gauge variables of Eq. (12. ip : the electron-photon, 
electron-ion and ion-photon rates of momentum exchange appearing in Eqs. (I2.8p - fl2.10p are 
given bj@: 



e2\2 8 /e2^2 



r^e = ^1-6(^67, r^i = 'T-iO'i^, ae7 = -7r( — j , o-i^ = -7r( — j , (2.11) 
Vei = ne\\ — (Tei = Tie, CTci = ;7;^lnAc, kQ = —\\- — . (2.12) 

V me i le* V ^eTT 

In Eq. (12. lip and (I2.12p . T and h are, respectively, physical temperatures and physical 
concentration^. 

2.3 The problem of the closure 

The two global variables entering the one-fluid description are the center of mass velocity of 
the electron ion system and the total current, i.e. 

m^Ve + rriiVi - ^ ^ . . 

Vh = \ , J = niVi-neVe. (2.13) 

rric + rrii 



^Note that Ac is the Coulomb logarithm PSII^ . 

''if the rates and the cross sections are expressed in terms of comoving temperatures T = aT and comoving 
concentrations n = n the corresponding rates will inherit a scale factor for each mass. For instance aFei 



becomes Uc \jT/{mea) (e /T ) In Ac, if comoving temperature and concentrations are used. 
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The one-fluid equations supplemented by the incompressible closure are often used as a 
starting point for the analysis especially in connection with the potential emergence of in- 
teresting scaling laws (see, e.g. [31]). The bulk velocity of the plasma coincides, in the latter 
case, with and the incompressibility condition dictates V ■ Wb = 0. Having said this, the 
reduction from the two-fluid to the one-fluid variables is essential, in the present context, 
exactly because V ■ "Ub 7^ and the velocity correlators depend on the density fluctuations 
of the plasma. The difference of Eqs. (12. 8 p and (12. 9 p leads to the evolution equation of the 
total current. Since the rate of Coulomb scattering is much larger than the conformal time 
derivatives of the current, the Ohm equation can be approximated by the standard form of 
the Ohm law (see fifth article of Ref. [H]): 

J = a(E + v^xB + ?^-:^^), (2.14) 

V eriQ erio J 

where uq = a^ftci denotes the common value of the (common) electron- ion concentration; a 
is the conductivity of the predecoupling plasma. After electron-positron annihilation cr(T) 
can be estimated as 



^(T) = ai — . ^-V, Ac(T) = Aflly^' = 1.105 xlO«^ 



0.02258 



(2.15) 

where cxi = 9/ {8n\/3) depends on the way multiple scattering is estimated. 

The sum of Eqs. (12. 8p and (12. 9p leads to the evolution equation of the bulk velocity of 
the plasma Vh which is directly coupled to the photon velocity and to the density contrasts 
of photons (i.e. 6^) and baryons (i.e. (5b): 



drv^ + nvy, = -W<i) + z/bWb + \^e\v^ - t/b), (2.16) 

a^Ph 3 Pb 



drV^ = -^V5^ - V0 + u^V% + e'({7b - v-y), (2.17) 

a.5b = Sdrtp - V ■ iTb + (2.18) 
4 - 

drS^ = Adr^ - -V ■ v^, (2.19) 

where z^b and z/^ denote the thermal diffusivity coefficients and e' = hcCr-yeO' is the differential 
optical depth of electron-photon scattering. Equations ( I2.16p - (l2.19p are all consistently 
written in the longitudinal coordinate system of Eq. (12. ip . 

Dropping the third and fourth terms in Eq. (I2.14p (i.e. the thermoelectric and Hall 
terms) the evolution equation of the magnetic fields reads, in the one-fluid approximation, 

drB = V X (4 X B) + XV^B, A = (2.20) 

Ana 
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If the Hall and the thermoelectric terms are kept in Eq. fl2.14l) . Eq. fl2.20p becomes 
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drB = V X (4 xB) + -^V'i? + V X (^] - -^V X [(V xB)x B]. (2.21) 



As in [17] we shall assume that the magnetic field is sufficiently small so that the Hall 
term is negligible (see, however, [M]). We shall also assume that there is no particular 
alignment between the gradient of the concentration and of the electron pressure so that the 
thermoelectric term is also negligible. 

The standard incompressible closure stipulates that the velocity field is solenoidal but 
this is not the situation described by Eqs. ( ]2.16p -( l2.20p where, in general, V ■ tly 7^ and 
V ■ 7^ 0. Prior to photon decoupling the electron-photon scattering rate drives the baryon 
and photon velocities to a common value which will be denoted by v^h (i-e. Wb — '^7 = v^h)- 
After summing up Eqs. fl2.16p and (12.171) to eliminate the scattering terms we arrive at the 
following equation 

where uth is the coefficient of thermal diffusivity of the baryon-photon system: 



i5(/?!+i)^^^' ^^'-h^' ^^-^^^ 



and Rhi^) is related to the sound speed of the plasma and it is defined as 

RM = -- = 0.629 ( -^^) ( ^\ c,b(r) = , ^ (2.24) 
^' 4p, U.02258ni09i; ^ J3[R^{t) + 1] ^ ' 



In Eq. (12.181) pb and 5b denote, respectively, the mass density of baryons and its inhomo- 
geneities: 

Pb = rnene + mini, 5b = (2.25) 

Pb 

In the ACDM case the vector modes of the geometry are absent. Strictly speaking the 
flow considered here is not only compressible (i.e. V ■ u^b 7^ 0) but also irrotational (i.e. 
V X v^b = 0). Since t7^b is irrotational it can be written as -iXyb = Vw^b- As an example the 
the equations for -u^b and the magnetic diffusivity equations can be written as: 

9..7b + ^^«7b = - ^^^^ - + -t.V u,b, (2.26) 

d^B = V X [Vn^b xB] + AV^fi, (2.27) 



^In the present analysis we shall assume a vanishing external gradient in the concentration of the charged 
species even if this requirement could be relaxed leading, presumably, to an effective evolution of the thermal 
and magnetic diffusivity coefficients with the physical scale (see, e.g. [26l [27]). 
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where, according to the notations and conventions of appendix |X] we used that: 



J X B = -a'^pj 
3 



(2.28) 



Equations fl2.26p and f l2.27p show that the incompressible closure is not consistent with the 
presence of scalar inhomogeneities; the adiabatic closure seems to be more appropriate and 
it is at least consistent with the hydromagnetic equations and with the ACDM paradigm. 
The adiabatic closure amounts to requiring that the fluctuation of the specific entropy (i.e. 
the ratio between the entropy density of the photons and the concentrations of the baryons) 
vanishes; this condition is equivalent to the requirement that 6^ = 4(5b/3. 

As in the case of the one-fluid hydromagnetic equations in flat spac^, different closures 
could be adopted (such as the isothermal closure or the closure with constant Ohmic current). 
While it could be interesting to explore other closures, the one suggested here is definitively 
better motivated from the viewpoint of the large-scale initial conditions imposed, at early 
times, on the Einstein- Boltzmann hierarchy [T5] . 

In the dynamo theory the compressible closure is usually adopted (see, e.g. [Il|2l[3]) and 
the same is true when the kinetic Reynolds number is very large (see, e.g. [32]). However, 
after electron-positron annihilation i?kin < 1 while Prmagn > -Rmagn 3> 1 • The hypothesis 
of primeval turbulence (implying the largeness of the kinetic Reynolds number) has been a 
recurrent theme since the first speculations on the origin of the light nuclear elements. The 
implications of turbulence for galaxy formation have been pointed out in the fifties by Von 
Weizsaker and Gamow [34] . They have been scrutinized in the sixties and early seventies by 
various authors [35] (see also [36] [37] and discussions therein). In the eighties it has been 
argued [3B] that first-order phase transitions in the early Universe, if present, can provide a 
source of kinetic turbulence and, hopefully, the possibility of inverse cascades which could 
lead to an enhancement of the correlation scale of a putative large-scale magnetic field, as 
discussed in [3ll [391 [10] • The extension of the viewpoint conveyed in the present analysis to 
earlier times (and larger temperatures) is not implausible but shall not be attempted here. 



3 Predecoupling flow 

The solution of the evolution equations fl2.16p - fl2.19p determines the correlation functions of 
the velocity field and the predecoupling flow. By taking the derivative of Eq. fl2.19p with 
respect to the conformal time coordinate r and by using Eq. fl2.22p a well known second- 
order equation for S.y{k, r) can be obtained and solved with different methods. In particular, 
within the WKB approximation [HI [121 [121 [B] , the solution for 6^{k, r) and u^h{k, r) reads, 

^In the context of flat-space magnetohydrodynamics the adiabatic closure amounts to a slightly different 
requirement, i.e. dr[pp^'^] = 0; pm is what we called pb (i-e. the mass density of the fluid) and k is the 
adiabatic index (i.e. the ratio of the heat capacity at constant pressure and constant volume) [29] . 
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in Fourier space: 

S^{k, t) = -^V'm(fc) + V^Mnik, r) cos [kr,{T)]e-'''/'''^ = Uh{k, r), (3.1) 

u^^{k,T) = -^AT7^(r)7^,(fc)sin[fcrs(r)]e-'^'/^^ (3.2) 

where rs(r) is the sound horizon and fcd(T) defines the time-dependent scale of diffusive 
damping, i.e. 

rsir) = r C3b(r') dr' = [ = I T 4(^0 dr' . (3.3) 

•^0 ^0 y3[i?b(r') + 1] 5 Jo 

The functions ip^nik^r) = T'ji{T)7l^{k) , M.-ji^k.r) and A^7e(r) 

-^7^(fc, ^) = ^ - 2)^^(r)7^.(A:), AT^(r) = -^=^^t^^T^(r), (3.4) 

are all defined in terms of T-ji{t) which is a simplified form of the transfer function of 
curvature perturbations discussed, for instance, in appendix B of Ref . [15] : 

TT^{r) = 1 - ^ r a\T')dT' = 1 - — /* a{t')dt'. (3.5) 
Jo a Jo 

Note that, for r — (where r^, denotes the last scattering) T7^(r^,) —(3/5). If A;r ^ 1 
we have, from Eq. (13. ip . that 5y — )■ —Sip^/S, Cgb and A^7e(r*) — )■ -\/3/5. Since the 

curvature perturbations are distributed as in Eqs. f l2.2p -( l23|) . the correlation function of 
the velocity for unequal times can be written as: 

{vM r) v,ip, t')) = ^U{q, \r - r'\) 5('\q + p), U{q, |r -r'\) = v{q) 5{r - r'), (3.6) 

where, to avoid confusions with vector indices, the subscript "7b" has been suppressed. The 
function v{q) appearing in Eq. (13. 6p is 

v{q)=r,V{q), V{q) =Ml{n)^-^Vn{q) sin^ [gr3(r,)] e"^^'/"!. (3.7) 

The correlation time Tc is, by definition, the smallest time-scale when compared with other 
characteristic times arising in the problem. Because of the exponential suppression of the 
velocity correlation function for r > Xd (where Td denotes the Silk time [30]), Tc approximately 
coincides with t^. The form of the correlator given in Eq. (13. 6 p is characteristic of Markovian 
conducting fluids [TTJ EI]- As a consequence of the smallness of Tc the velocity correlator for 
unequal times can be approximated with its Markovian form given in Eq. (13. 7p . The velocity 
field is exponentially suppressed for r > Tc = ^ / {k'^.^-^Vth) where Tc ~ 1.95x 10~^(z* + l)^ Mpc, 
z^. ~ 1090 and, typically, fcmax — ^d- Note that v{q) contains Tc and it does not have the 



9 



same dimensions of V{q). Defining, for immediate convenience, V^iq) = v{q)/{2TT'^), the 
velocity correlator can be expressed in physical space as 



{vi{x, t) Vj{y, t')) = \ VT{r)6ij + 



li(r)-VT(r) 



S{t-t') 



Vr(r) 



dq 



q3^3 



sm qr — — ■ cos qr 



VL(r) = / ^n(g) 

J q 



2 2 sin qr 
— - sin qr + — -— cos qr H 



(3.8) 



(3.9) 



where r = \x — y\. Thanks to the explicit form of the fluid flow, Eqs. f l2.2ip and fl2.27p can 
be solved by neglecting all the terms which are of higher order in the magnetic field intensity 
and by focusing the attention on the coupling of the compressible flow to the magnetic field; 
by following Ref. [HI EI] Eqs. ^I7n) - ^m\ can be solved iteratively as 



Bi{k,T) 



B. 



(n+l) 



(fc,r) 



£ Gkir - n) dn J d'q J d'p 6^'\k - q- p) 



n=0 

(-^ 



(27r)3/2 Jq 

X emni (^abn {qm + Pm) Va{q, n) B^f^\p, Ti) 



(3.10) 



(3.11) 



where, for simplicity, A is assumed to be constant in time; moreover, as in Eq. (I2.27p . 
Va{x, t) = dau{x, t) . The magnetic field can then be averaged over the flow by using either 
Eq. (13.61) or Eq. (13. 8p (depending if we work in real space or in Fourier space). From Eq. 
(13. lip the first few terms of the recursion are B\^\k^T), Bl^\k,T) and Bl'^\k,T). The first 
term is Bf'\k,T) = Qk{T)Bi{k) where Bi(k) parametrizes the initial stochastic magnetic 
field. Denoting with Hi{k,T) the magnetic field averaged over the fluid flow, the terms 
containing an odd number of velocities will be zero while the correlators containing an even 
number of velocities do not vanish i.e. (sf^"^"*^^) = Hf''^^^^ = and = uf'^^'^^ ^ 0. 



So, for instance, (sf ^) = while (5f ^) = iff ^ is 



(2)\ 



r(2) 



H?\k,r) 



(2vr)= 



d'^q / d^p / d^q' / d^p' 6^^\k-q-p) 6^^\p-q' -p') 



dnGkir - Ti) / dT2Qp{Ti - T2) {qm +Pm) {q'm' +Pm')^hm'n' b' n' emnieab, 



X {Va'{q\ T2) Va{q, Ti)) By{p). 



(3.12) 



After averaging the whole series of Eq. (I3.10p term by term the obtained result can be 
resummed and written as: 



H,{k,T) = {Bf\k,T)) + {Br\KT)) + {Br\k,T)) + ... = e-'''-B,{k) 



>(2) 



(3.13) 



where the magnetic diffusivity coefficient A = 1/ (47rcr) inherits a modification stemming from 
the large-scale flow: 



A = A + fo. 



-:—2 , , t dk 



vo = -iMj^in) / -^Vn{k) sin' {k/k,)e 



2 /7„2 
d 



(3.14) 
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and = l/rs(r^,). The averaging suggested here has been explored long ago in the related 
context of acoustic turbulence by Vainshtein [21]. Prior to decoupling, however, both (t>^) oc 
A-ji <^ 1 and i?kin ^ 1- The iterative solution indicated in Eq. fl3.13p seems then to be 
better defined, in the present case, than in the a kinetically turbulent plasma with strong 
inhomogeneities. The presence of Vq defines a new diffusion scale associated with the large- 
scale flow and this scales can be estimated from Eq. fl3.14p as 




where subscript M stands for "Markovian" and where denotes the (comoving) angular 
diameter distance to last scattering: since ku depends on time, it has been evaluated at last 
scattering. In the appendix [B] the analysis in the non-Markovian approach will be outlined, 
under some simplifying approximations, with particular attention to the determination of 
the diffusive scale. 



4 Evolution of the power spectra 

Instead of averaging the magnetic fleld intensities, as explored in the previous section, it is 
more practical to study directly the evolution equations of the magnetic power spectra with 
the aim of connecting them to the two-point function of the velocity fleld. In Fourier space 
the magnetic power spectrum can be deflned as 

{B,{k,T)Bj{k',T)) = Mi,{k,T)6'^^\k + k'), M,,{k,T) = P,^{k)M{k,T), (4.1) 

where Pij{k) = {5ij — kikj). The corresponding power spectrum in physical space is given by 

M,,(r, r) = MT(r, r)5,fc + [Mi^{r, r) - M^ir, r)]^, (4.2) 

^ + -(Ml-Mt) = 0, (4.3) 
or r 

where Eq. (14.31) is a consequence of the transversality condition, i.e. diM^^ = 0. The 
evolution equation of the magnetic power spectrum in the compressible limit stems from Eq. 
fl2.27p and can be obtained by taking the conformal time derivative of Mij{k,T) as deflned 
in Eq. (14. ip . The result of this step is the sum of two terms: each of the terms contains the 
product of the time derivative of the magnetic fleld with the magnetic fleld itself for different 
wavenumbers. The time derivative of the magnetic fleld can be eliminated by means of Eq. 
(I2.27P while the magnetic fleld itself can be written in terms of the integral equation derived 
always from Eq. (12.271) . The result of this procedure is, in the compressible case, 

/d'^Q r d'^cf' 
^2^^3/2 J ^27]-^3/2 ^rnkm' ^imn^abn^jm'n'^a'b'n' ^ 

X {vM T)vaiq,T')) M,b'{p,p, r, r '), (4.4) 
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where p = {k + q) and p' = {k + q'). Using now the exphcit form of the velocity correlator 
given in Eq. (13. 6 p and taking the trace of both sides with respect to the free tensor indices 
i and j the following integrodifferential equation is readily obtained 



drM + 2k'^\M 



d^q 
(2^ 



v{q) M{p,t) 



(k-p)^ Jk ■ q){k ■ p){q ■ p) 



p^ 



+ 2 



(4.5) 



where, for convenience, p has been kept explicitly but it is in fact p= {k — q). Note that A 
is the diffusivity coefficient accounting also for the effects of the flow (see Eq. (13.141) ). 

The same derivation can be performed in the incompressible case but bearing in mind that 
the correlation function of the velocity field will have a different form which is determined 
by the transversality of the fiow. The analog of Eq. (14. 5 p in the case of the incompressible 
closure is given by: 



drM + 2k'^\M 



cP'q 
(2^ 



v{q) M{p,t) 



{k ■ q){k ■ p){p ■ q) 



(4.6) 



which coincides with the equation derived by Kazantsev [T7] except for the presence of 
l/(27r)'^ coming from the different conventions on the Fourier transforms (we follow here 
the conventions spelled out explicitly in appendix |A|) . To correctly interpret Eq. (14. 6p it is 
essential to point out that the velocity field parametrized by v{q) is now solenoidal (rather 
than irrotational as in the case of Eq. (I4.5P ). To derive Eq. (14. 6 p we used the analog of Eq. 
(13. 6 p but valid for the standard incompressible closure, i.e. 



(fi(g,r) Vj{p,T')) 



5r 



U{q,\r-T'\)5^'\q + p), U{q,\T - t'\) = i{q) Sir - r'). 

(4.7) 



The results of Eqs. (14. 6 p and (14. 7p are a consistency check of the whole approach. For 
sake of comparison, the incompressible equations (derived from Eq. (14. 6p ) will be sometimes 
reported after the corresponding results valid in the compressible case. It is finally useful to 
point out that in Eqs. (14. 5 p and (14. 6 p the diffusivity coefficient is A (and not A itself). This 
property has been already derived, at the level of the field intensities, in Eqs. (I3.13p -( l3.14p . 
Strictly speaking Eqs. (14. 5 p and (14. 6 p describe the evolution of the magnetic power spectra 
averaged over the fiuid fiow. In this respect following the notations of Kazantsev [T^ (see 
also |2T]) it is sometimes useful to distinguish the "double" stochastic average (( ... )) (over 
the velocity and over the magnetic field) from the single average (valid either for the velocity 
or for the magnetic field at the level of the correlators). To preserve a certain simplicity in 
the notations the double stochastic average has been avoided. 



4.1 Diffusive approximation 

The integrodifferential equations (14. 5 p and (14. 6 p can be studied, with two complementary 
approaches, either in Fourier space or in physical space. By taking the limit g — )■ (while k 



12 



is held fixed), Eqs. fl4.5p - fl4.6p lead to a diffusion equation which is very similar to the one 
often encountered is the dynamics of wave turbulence [HIIIT]. Equations (I4.5p - fl4.6p can 
also be transformed into a Schrodinger-like equation in physical space where the analog of 
the wavefunction is related to the power spectrum introduced in Eq. (14. 2p . In what follows 
the relevant evolution equations will be derived in these two complementary approaches. 
The attention will be focused on the compressible mode, however, as a cross-check, the same 
derivations have been also performed with the incompressible closure. In the latter case 
the equations reproduce exactly the results reported in [T7]. This represents an important 
cross-check for the consistency of the whole procedure. 

In the diffusive limit Eq. (14. 5 p is expanded for e = \q/k\ -C 1. The integrand appearing 
at the right hand side of Eq. (14. 5 p written in terms of e becomes 



M(p, r) 



2 {k.p)^ 



ik ■ q){k ■ p){q ■ p) 



k'^ui kVl - 2ex + 



q2p2 

(1 - exy 
1 - 2ex + e2 



x(l — ex){x — e) 



1 - 2ex + e2 



(4.8) 



where x = cos^ and ^ = {q-k). After expanding Eq. (14.81) in powers of e, direct integration 
over X between —1 and 1 gives 



3 



2A;2r BM .d'^M 

2M + 6k-— + 3k^ 



15 



dk 



dk^ 



(4.9) 



Using the resuhs of Eqs. dMD and (gS]), Eq. (gS]) 



dM 
'd7 



+ 2\k'^M = 7 



2M + 6k— + 3k^— 



V2 



d?q 
3\J (2^ 



(4.10) 



where 7 = 2f2/5. From the explicit expression of v{q) it is possible to estimate 7 by replacing 
sin^ -T- 1/2, A^7^(t^,) — 7- 3/25, and by fixing the upper limit of integration at kd_: 

An 



^ = 30kTT)^ [k-J ^p"- 



V = IT, 



(4.11) 



where, as already mentioned, Ai (r*) — )■ 3/25; the introduction of the variable rj will be 
useful in section |5] since it will simplify the expressions of the Laplace transforms. Choosing, 
for sake of simplicity fed ~ 0.1 Mpc~^ as well as the fiducial set of parameters of Eqs. (12. 3p 
and (12. 4p rj can be explicitly estimated; for r ~ 0{tj,) we have that rj ~ 6.04 x 10~^^. 

As anticipated after Eqs. (14. 60 - 04. 7p . the same kind of equation obtained in the com- 
pressible limit can be derived when the fluid flow is incompressible. Indeed, from the right 
hand side of Eq. ( 14.60 we have 



M(p,r) 



k' 



{k ■ q){k ■ p){p ■ q) 



k'^Mi k^l - 2ex + 



x{l — ex){x — e) 



1 - 2ex + 



• (4.12) 
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By expanding Eq. fl4.12l) for small e and by integrating over x between —1 and 1 we have 



4,2.. 2P 

-k^M+ 

3 15 



4M + 2A; 



dM 
~dk 



implying, from Eq. (14. 6p . 

-— + 2\k^M = -V2 
or 5 



AM + 2k-— + k^ 



dk 



dk^ 



dk^ 



^2 = TTT 



3! 7 (27r 



(4.13) 



^ ^ g^S(g). (4.14) 



Apart from the slightly different conventions, Eq. (I4.14p reproduces the corresponding results 
of [T7] and corroborates the correctness of the results obtained here in the compressible case. 



4.2 Schrodinger-like equations 



In physical space Eq. (14. 5 p can be transformed into a Schrodinger-like equation whose 
generalized wavefunction is related to the magnetic power spectrum. Since, by definition. 



M(r,r) 



(27r)3/2 



d^ke-'''-''M{k,T), 



(4.15) 



both sides of Eq. (14. 5 p can be multiplied by exp [—ik ■ r] and then integrated over d^k. 
Recalling that k = q + p, after simple algebra the evolution of M(r, r) can be obtained in 
physical space: 



drM - 2A 



dradr°- 





-2 




GijV 




GikFj^j 



(4.16) 



where it is understood that M = M{r, r). According to Eq. (13.80 . v coincides with the trace 
of the correlation function of the velocity in real physical space, i.e. v = (2Vt + ^l)- The 



tensors Gij and Fij are: 



Gij = M5,j + (M - 3M) 



r2 ' 2 
where M(r, r), M(r, r) and ML(r, r) are related as follows: 



Sij + (3Vl - V) 



2r2 



M(r, r) 



Since for a generic function /(r, r) it is easy to show that 



1 d 



dvadr"- dr 
Eq. ( ]4.16p can be transformed as: 



ML(r,r) 1 r ^ 
= ^ / X M[x,t) ax. 



1 



dr . 



dvidrj 



Virjf] 



V]. 



dM -B'^M 
+ 2(\/L-A)^r^+2 



dr 



or r 



dM 



dr 



+ 2 



(4.17) 



(4.18) 



(4.19) 



M = 0, (4.20) 
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which can also be written in the following equivalent form: 



dM 2d 
Or Or 



r\V^ - A) 



dM 



dr 



+ 



or r 



M = 0. 



(4.21) 



By performing the same derivation in the case of Eq. f l4.6p . the analog of Eq. (14.2 ip becomes 

(4.22) 



dM -^d'^M , _ _15M 2—dv 

-- + 2 V - A — - + 2[v - 4A + V)-^ + -M— = 0, 
or or-^ r or r or 



where, in this casqj 



V = ^ I x\{x) dx = ^, 2Vt + Vl = v. 



r^ Jo 



(4.23) 



Equation f l4.23p coincides with the analog equation obtained in [T7j . It is finally appropriate 
to mention that Eq. fl4.2ip coincides with the equation derived by Vainshtein and Kichatinov 
in [IH] once the relevant functions are appropriately renamed. The correspondence is such 
that Vt Tnn, Vl ^ll, and 5ll Ml = 2M. Equations fOTj) and dM]) can be put in 
different but equivalent forms like, for instance, 

dM 

+ 2£;(r)M(r,r) = 0, 



dM 



+ 24(r)M(r,r) =0, 



(4.24) 
(4.25) 



where Eq. fl4.24p holds for an incompressible flow while Eq. f l4.25p holds for a compressible 
flow. In the incompressible case the linear operator Cf{r) takes the form 



^fir) = -f{r] 



dr"^ 



/i(r)|: + /2(r). 



where the three functions /(r), /i(r) and f2{r) are defined as 

fir) = A - F, A(r) = i[4A - F - /,(r) = 



Conversely, in the compresible case, the linear operator jC,g{r) is 



d^ 



d 



and the explicit form of the functions g{r), gi{r) and g2{r) is 

9{r) = A - VUr), gi{r) = |-(A - li.) + ^(A - li), 

92{r) = -^{v + \/l) + 4(3^L -v) = -^{Vt + Vi) + ^{V^ - Vt) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



^We recall that, in the incompressible case, the explicit forms of the velocity correlator in real space is 
the same as the one given in Eq. p.Sp but with the difference that the transversality condition applied to 
the velocity field implies that drVh + 2(Vl — Vt)/?" = 0; the latter condition leads to the first expression of 
Eq. (14:2311 . 
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By focusing the attention on the compressible case of Eqs. fl4.25p and (I4.28p - fl4.29p . Eq. 
fl4.25p can be rewritten by ehminating the first derivative with respect to r: 

d^^ — 

-2iy(r)^, M{r,T) =Q{r)'${r,T)] (4.30) 



9^ _ , 



the potential W{r) and Q{r) are 



W{r) = g 



2\g 



A\g 



f ^ Q' 9l 



(4.31) 



where the prime denotes, in Eqs. f l4.3ip and fl4.32p a derivation with respect to r. Defining, 
for immediate convenience, gi/g = F'/F, Eq. f l4.3ip implies, together with the definition of 
^i(r) of Eq. flCT]) . 



W{r)=g{ry-^^ + g2{r), g(r) 



F{r) 



F{r) = r^^). 



(4.32) 



5 Solutions for the compressible mode 



5.1 Kinetic equations for the power spectra 

It is interesting to solve the equations obtained in the previous section by assuming that the 
magnetic power spectra are assigned at a specific pivot scale conventionally denoted by /cl 

m 

27r2 



-VBik, T)Pij{q), Pij{q) = 6ij - qiqj, 



(5.1) 



where qi = qi/q. The evolution equation of Pb(^, t) will be solved by imposing the following 
boundary conditions 



P(A;,0) = Ab(-' 



V{K,t) = AbG{t). 



(5.2) 



Even if G{t) will be left unspecified a simple choice would imply G(r) cx 5(r). The evolution 
equation for V-^ik^T) stems from Eq. fl4.9p by comparing the parametrization of Eq. (14. ip 
with the one of Eq. (15. ip : 



dr 



20 VB-12k 



dk 



+ 3k' 



dk^' 



(5.3) 



Depending on the value of the conductivity there are two regimes which can be identified. 
In the first regime the magnetic diffusivity is subleading (i.e. k'^X -C IO7) and Eq. (15. 3p 
becomes: 

&Pb 



dr 



7 



20 VB-I2k 



dVB , _ 2 d^'PB 



dk 



+ 3k' 



dk^ 



(5.4) 
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When the magnetic diffusivity is subleading Eq. f l5.4p can be solved with the boundary 
conditions fl5.2p by standard Laplace transform methods. In particular, defining 



(5.5) 



Eq. (15. 4p together with the first of the two conditions of Eq. (15. 2p implies the following 
(ordinary) differential equation: 



+ 



20 

3P 



A, 



3'jk'^ \ kL 



k \"B-i 



dk"^ k dk 

After a standard change of variables Eq. (15. 6p can be transformed as 

1 I 5~ 

37 



d^Qi 



g2(.,7)QB = -4^e^^ 



where the variables z and Q-q{z^s) are simply: 

z = \n (k/ki^), Q^{z, s) = e-^^/2pB(^, s). 
The general solution of Eq. (15. 7p is 



p = nB- -. 



(5.6) 



(5.7) 



(5.8) 



Qb{z,s) 



ciis) + 



At 



A, 



67g(p- g(s,7)) 



6'^' + 



A, 



C2(S - 



67g(s,7)(p + g(s,7)) 



-g(s,7)z 



37 — g^(s, 7) 
If 2; > the second condition of Eq. (15.20 implies that 

1 



(5.9) 



C2[S) 



Ar 



+ 



1 



67g(s,7)(p + g(s,7)) 37(^2-^2(5,7)) 
Ab 



Cl[S) = r-, (5.10) 

67g(s,7)(p-g(s,7)) 

where g{s) simply denotes the Laplace transform of the function G{t) which appears in the 
second relation of Eq. (15. 2p . Equation (I5.10p applies for z > 0, i.e. k > ki^. Similarly, if 
z < 0, the second condition of Eq. (15. 2p together with the requirement of the existence of 
the inverse Laplace transform implies 

1 1 



A, 



9[s) 



67g(s,7)(p- g(s,7)) 37(^2 - q'^{s,-f)) 
Ab 



C2(S) = 7 77 7 TT. (5-11) 

67g(s,7)(p + g(s,7)) 

Equation (15. lip applies for ^ < 0, i.e. for k < ki^. The Laplace antitransform can be easily 
deduced by appropriately choosing the integration contour in the complex plane. In the case 
2 < we have that Qb{.z, s) is given by: 



Qb{z,s) 



Ab 
37 



.p2 _ q^{s,'^) 



+ 375((s) e' 



9(^.7)2 _ 



p2 _ g2(s^^)_ 



(5.12) 
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An analog expression can be deduced in the case z > by choosing the constants of Eq. fl5.9p 
as in Eq. f lS.lOp . Since the exphcit expressions of g(s,7) and p have been given in Eq. (15. 7p 
it is easy to see that the first term in Eq. (I5.12p has a branch point from the argument of 
the exponential and a simple pole from the denominator; the second term has only a branch 
point while the third term has just a simple pole. By taking the inverse Laplace transform 
of Eq. ( I5.12P the spectrum can be explicitly obtained. The same procedure described so far 
can be repeated in the case z >0. The general form of Vb{z,t) is then given by: 



A, 



+ 



\z\e 



5z/2 



y/l2n^u 



^5 7m/4-|z|2/(127«) 



du 



u 







1 






Jo 



sm 



37 ) 



■ dx. 



(5.13) 



(5.14) 



x + 37(nB- 7/2)2 

The first and second terms of Eq. f l5.13p come from the contributions containing, respectively, 
a simple pole and a branch point in the Laplace transform Qb{Zi s)\ the third term in Eq. 
fl5.13p corresponds to the piece containing a simple pole and a branch point. The properties 
of the solution fl5.13p are determined by the value of 7r already examined in Eq. (14. lip . It 
is useful to note that when 7r is kept fixed in such a way that 7 — )■ and r — 00 Eq. (15.140 
has a definite limit: 

(5.15) 



lim J-'(r, ne) = Erf 

7— >0, r— >oo 



,v/T2^. 

As already stressed in connection with Eqs. (I5.3p -f l5^ . the solution (I5.13P holds in the 
limit k^X < IO7. In the opposite limit (i.e. fc^A > IO7) it is possible to look for the solution 
by bootstrapping Eq. (15.130 : 



VB{k,T) = e5 7r/4+3(nB-7/2)2 7r g^f^y 



(5.16) 



Inserting Eq. (15.160 into Eq. (15.40 and eliminating the first derivatives with respect to k the 
following equation can be obtained in terms of the rescaled function = Q{k,)/k'^: 



dk"' 



+ 



2A 

"3^ 



1/4 



fc2 



0, 



z/ = |nB-7/2|. 



(5.17) 



The general solution of Eq. (15.170 is given in terms of a linear combination of modified Bessel 
functions Iu{z) and Ky^v) [49j. By assuming that the initial magnetic power spectrum is 
such that ub < 7/2 the solution of Eq. (15.170 is: 



G{k) 



Ar 



(5.18) 



which satisfies the correct boundary conditions since, in the limit k ^ k^, we have that 
G{k) — )■ AB^k / ki^)"'^~^ . The results obtained so far partially challenges the standard argu- 
ment leading to the calculation of the magnetic diffusivity scale which is phenomenologically 
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important insofar as it determines the scale of the exponential suppression of the magnetic 
power spectrum. Prior to decoupling the magnetic diffusivity scale is estimated by requiring 
that k"^ ~ 47r(T'H. The latter condition stems directly from the magnetic diffusivity equation 
and it totally neglects the flow. 

The solution (15.131) can be used to compute the evolution of the spectral index induced 
by the predecoupling flow. The spectral index of the magnetic field is defined as 

Having assigned the spectrum at a fiducial (pivot) scale, it is interesting to ask what happens 
to the power spectrum for larger length-scales (i.e. smaller wavenumbers) . From Eq. f l5.13p 
it is possible to derive an evolution equation for as it is defined in Eq. (I5.19p . After simple 
algebra the result is given by: 



dN 
'd7 



7 



-15- 



dN 
dz 



/dN\ 



dSN 



The solution of Eq. fl5.20p can be found in the form N{z, t 
conditions 6N{z, 0) = 0. In this case Eq. fl5.20p becomes: 

ld6N _^d^6N ^^^^ 
7 dr dz"^ dz 

Using the standard Laplace transform technique we have that 

, . 1 



(5.20) 

riB — 1 + SN{z, t) with initial 



N{z, r) = TiB + 



z + 67r ub — 



7 



(5.21) 



(5.22) 



127r73r3 ' '''''' ^/W^ 

The result of Eq. (15.220 matches pretty well with the result obtained directly from Eq. 
fl5.19p by using the exact solution of Eq. (I5.13p . as it can be seen by plotting the respective 
functions for specific values of the spectral index. To avoid lengthy digression this analysis 
will be omitted. 



5.2 Evolution in physical space 

The Schrodinger-like equation in the compressible case (see Eqs. (I4.30p . ( I4.3ip and (I4.32p ) 
allows for a semi-quantitative description of the limits r (where magnetic diffusivity 
dominates) and r — > oo (where the large-scale flow dominates). To deduce the potential 
in the limit r — i- the correlation functions of the velocity can be expanded in the limit 
R = rk^ < 1 where k^ is an auxiliary scale in the space of the wavenumbers. Let us therefore 
start with Vt(?") and VL{r) written as 

V^{r) = MlirMnC^y T —x''-'sm^xa)e-'^"^A{x,R), 

\ Kp J J ko/kp X 

li(r) = Mnin)An(j^) / — x^-^in^ (x«)e-2-'^ 5(a:, i?), (5.23) 
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where /cq can be estimated from the inverse of the comoving angular diameter distance to 
last scattering. In Eq. fl5.23p the functions A{x,R) and B{x,R) are defined as 



A{x,R) 
B{x,R) 



sin xR — -- cos xR 



X- 



'R^ 



2^2 „ smxR 

sm xR H — ^^-77 cos xR + 



x^R^ x2i?2 ' xR 

In Eq. fl5.23p the following rescaled quantities have been introduced: 

The expansion of Vt(-R) and Vl(-R) for i? < 1 becomes therefore 



R = rk^. 



where 



VT{R) = Cin,,An) 
Vj^{R)=C{n,,An) 

C{ns,An) = T~c—{ -r 



R^ R^ 

XK, a, 13) - XK + 2, a, /?)— + X{n, + 4, «,/?)— + ... 

X(n„ a, /3) - T{n, + 2, a, + X(ns + 4, a, + ... 

lU 00 



25 Vfe 



n-l 



X{ns,a,l3) = — X 
Jo X 



sin^ (xa)e-'^"f'. 



(5.24) 
(5.25) 

(5.26) 
(5.27) 

(5.28) 



To discuss in a semi-quantitative manner the properties of Eq. fl4.30p it is useful (even if not 
necessary) to identify ~ n Idxij^). Following the approach of Zeldovich (see, e.g. [3]) it 
is useful to introduce an interpolating form of the correlation function of the velocity which 
is suppressed for r — )■ cxd and leads to the correct limits of Eqs. f l5.26p and f l5.27p for r — )■ 0. 
The interpolating form is given by 



VAr) 



voe 



(5.29) 



where Vq has been already introduced in Eq. fl3.14p and coincides with the normalization 
which can be derived from Eqs. fl5.26p - fl5.28p in the limit Ai-jzi'T*) ~^ \/S/5. Using Eq. 
fl5.29p the effective potential of Eqs. f l4.3ip - fl4.32p can be computed. For instance g{r) = 
vo[l + 6 — exp (— 3yur2)] where e = X/vq <^ 1. In the limits r — )■ and r — t- oo the potential 
W{r) becomes 

2A 2A 
lim W(r) = — , lim W(r) = — , (5.30) 

showing that for intermediate and large scales the suppression of the correlation function of 
the magnetic field is determined by the large-scale flow contained in A. If f o ^ A the potential 
can become negative. Equation fl4.30p may have negative "energy" levels, i.e. solutions for 
\E'(r, r) ~ exp [— 2£'r]$(r) which are exponentially increasing in time (i.e. E < 0). For at 
least one level to exist we should have, roughly, that 



ri 



r.^E-W{r) n 
dr>-, 



gir) 



(5.31) 
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where ri and r2 are the inversion points. Equation fl5.3ip can be verified but the growth rate 
of the field is F ~ fivo leads to a neghgible integrated growth at last scattering, i.e. Tr^: ~ 
C(10~^^) even assuming y/JI ~ l/dA^T^)- This semi-quantitative argument (also employed in 
the case of gyrotropic turbulence [3]) follows from the analogy with the Schrodinger equation 
in the WKB approximation. 

6 Concluding remarks 

The standard treatment adopted for the evolution of predecoupling magnetic fields usually 
neglects two aspects which are the starting point of the present investigation, namely the 
compressibility of the plasma and the large-scale fiow induced by curvature perturbations. 
A third related coincidence is that, prior to decoupling, the magnetic Reynolds number is 
roughly 20 orders of magnitude larger than the kinetic Reynolds number which is, in turn, 
smaller than one. The basic idea has been to derive an integrodifferential equation valid for 
the compressible mode and for a standard adiabatic closure. In Fourier space, the diffusive 
approximation leads to a nonlocal diffusion equation similar to the ones often discussed in 
wave turbulence. In physical space the integrodifferential equation can be transformed into 
a Schrodinger-like equation whose effective potential depends on the spectral properties of 
the large-scale flow. Some applications have been discussed by solving the corresponding 
equations: they range from the effective evolution of the magnetic spectral index to the 
qualitative discussion of the spectrum of the eigenvalues of the equation for the magnetic 
power spectra in physical space. 

In the absence of large-scale flow the magnetic power spectra are suppressed at small 
scales (i.e. large wavenumbers) because of the flnite value of the conductivity which implies 
an effective (comoving) wavenumber ka 



where — 14116 Mpc denotes the (comoving) angular diameter distance to last scattering. 
In units of we have that k„dA — O{10^^) implying that the corresponding length-scale 
is much smaller than the Hubble radius at last scattering, as expected. The results of 
the present paper suggest, however, that the correct magnetic diffusivity length-scale is 
much larger than k~^. Using the Markovian approximation it has been shown that the 
magnetic fleld can be averaged over the large-scale flow and the resulting diffusive scale can 
be estimated from Eq. (13.151) and it is 



where the flducial set of parameters of Eq. (12.41) has been assumed in the context of the 
vanilla ACDM scenario. Not only /cm ^ k„ but, in units of d\, it turns out that /cmc^a — 




(6.1) 




(6.2) 
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(9(10^)(iA- From the analysis of the evolution equation of the power spectrum in the diffusive 
approximation it is possible to derive yet another dissipative scale, i.e. (see Eq. flS.lSp 
for a definition) whose explicit value, always in terms of our fiducial set of parameters, is: 



By comparing Eqs. (16. ip . (16. 2p and (16. 3p it is clear that the following (approximate) hierarchy 
holds, i.e. ky[ < k^ < k„. Both the results for /cm and k^ have been derived, directly or 
indirectly, by assuming the Markovian approximation for the velocity correlator. If the 
Markovian approximation is relaxed the diffusive wavenumber does not get larger but even 
smaller. From Eq. flB.22p it is possible to obtain that 



In summary, it has been shown that the large-scale flow can affect the evolution of the mag- 
netic power spectra at large scales not only by potentially shifting the effective spectral index 
but also by changing the diffusive scales. Thanks to the results derived in the present paper, 
the evolution of the magnetic power spectra prior to decoupling can be addressed in terms 
of a novel set of equations which take into account the effects of the large-scale flow and 
which can be explicitly solved in various physical limits. According to the present results it 
does not seem correct to treat the predecoupling plasma by simply assuming an incompress- 
ible flow; the latter closure is sound in the absence of large-scale curvature perturbations, 
in fiat space-time and for high (kinetic and magnetic) Reynolds numbers. None of these 
three assumptions are verified after electron-positron annihilation and prior to last scatter- 
ing. A closer scrutiny of the description developed here seems therefore both motivated and 
potentially rewarding. 
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A Basic conventions 



In this appendix some basic conventions will be summarized. As an example, the Fourier 
transform of the magnetic field and of the velocity field are defined as 



Bi{x,T) 



cP'q 



(27r)3/2 



(27r)3/2 



(A.l) 



While the magnetic field B is strictly divergenceless (i.e. V ■ -B = 0), the velocity field may 
be either solenoidal (i.e. V ■ "U = 0) or not solenoidal (i.e. V ■ v ^ 0). If the flow is strictly 
solenoidal it is also incompressible. In the vanilla ACDM case the flow is irrotational and 
compressible. The following definitions will also be employed in section |21 



QB{k, t) e 



-ik-x 



(27r)3/2 

where J1b(^, t) and liijik, r) are given by: 



Ilij{x,T) 



d^k 
(27r)3/2 



Ilij{k,T) e 



-ik-x 



fiB(fc,r) 
Uij{k,T) 



d^q 



TTa^p^J (27r)3/2 



Bk{q,T)B\k~q,T), 



d^q 



47ra4 J (27r)3/2 [ 



Bi{q, T)Bj{k - q, r) - -Bk{q, T)B\k - q, t)S, 



Recalling that (p^ + p^)V^(Tb = didjU^Q we also have, in Fourier space, that 



-> 3 ki kj 

(jB(fc,r)- 



d^q 



Wna^p^J (27r)3/2 
Finally Eqs. flA.2p and (lA.SP imply 



B\q, T)B^{k - g, r) - —Bk{q, T)B\k - g, r) 



B\q,T)B^{k-q,T). 



(A.2) 



(A.3) 
(A.4) 



(A.5) 



(A.6) 



B Non Markovian approach 

In the non-Markovian approximation the correlator of the velocities in the compressible case 
is still Gaussian and it is proportional to the correlator of the curvature perturbations: 

{v.iq,n)v,{p,T2)) = ^T{q,ruT2){nMn,{p)), (B.l) 

having introduced the function T(q, Ti, T2) 

T{q,Ti,T2) =Mn{ri)Mn{r2) sin(gc,bri) sin (gcsbr2) e-''''^"^(^^+^^\ (B.2) 
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where A^7^(r) is given by 

^ ^ y3(l + i?b)3/^ ^ 5 ^ ^ 

The hmit in the second expression holds after matter-radiation equahty and for i?b 0. As 
explained this approximation is justified at last scattering where i?b ~ 0.6. The expression 
of Eq. (]B.2[) can also be written as: 

r(g, n, ra) = ^jcos [qc^^in - r^)] - cos [qc^^in + r2)]}e-'?'^«^(^i+^^)7W7e(ri) Mn{r2). (B.4) 
With these notations the correlator of the velocities can be written as: 



{vMri)vMr2)) = ^viq)Tiq,n,T,)6^'\q + p), v{q) = %Vn{q). (B.5) 

To simplify the problem it is practical to use the limit a — ?■ oo where the magnetic diffusivity 
equation is given by: 

_ = V X (tJ^b X B). (B.6) 
By writing Eq. (1B.6P in Fourier space we have: 

drBi = ^^^^1/2 J J d^P^'-^Kk - q- P) ernni^abniqm + Pm)Va{q,T) Bb{p,T), (B.7) 

The solution of Eq. (]B.7|) can be formally written as: 



Bi{k,T) = ^3^/2 J j d^P^^^Kk- q-p)emni(^abn{qm+Pm)Va{q,Ti)Bi,{p,Ti). 

(B.8) 

Equation (iRSj) can be solved by iteration as 



oo 



n=0 



Bt^'\k,r) = ^i-^ £ dn J d\ J £p 5^'\k-q-v) 

X emnieabn{qm+Pm)Va{q,Ti) BI'^\p,Ti). (B.9) 

We can then average the magnetic field over the velocity field. The terms containing an 
odd number of velocities will be zero while the correlators containing an even number of 
velocities do not vanish, in formulas 

^^(2.+i)^ = i/f = 0, (i?f = iff ^ 0. (B.IO) 
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This conclusion holds both in the Markovian and in the non-Markovian case since it is directly 
related to the Gaussianity of the curvature perturbations. Let us now write explicitly the 
contributions up to second order and let us compute the first few terms in the expansion: 



d'q' / dy S^'\p-q'-p') 



(27r)3/2 

X ebm' n' Ca' b' n' {q'm' + P'm') I ^'^2 Va' T2) Bb{p), 



(B.ll) 



(B.12) 



where we used that B^\p',T2) = Bb{p'). Inserting Eq. (1B.12P inside Eq. (IB. lip we have 
that 



Bl'\k,T) 



(2vr)= 



d^q / d^p / d^q' / dy 5^^\k-q-p) S^^\p-q' -p') 



X 



/ dri / dT2 {qm + Pm) iqm' + Pm')^b 

Jo Jo 



X Va' {q , T2) Vaiq, Ti) By (f) . 

By averaging over the velocity fields we have that: 

H^'\k,T) = 1^ Jd'q j d^p j d\' j d^'p' 5^^\k-q-p) 5^^\p-q' -p' 



(B.13) 



X / dri I dT2{qm+ Prn) iq'm' ^ Pm')^bm'n' ^a'b'n' ^mni^abn 

Jo Jo 

X {va'iq', T2) Vaiq, n)) By{p'), (B.14) 

where {Bl (/c,r)) = Hi {k,T). After inserting the explicit expression of the correlator 
obtained in Eq. flB.51) . (Eq. flB.MD reduces to 

nPik, r) = J^Jd'qJ d'p J d\' J d'p' 6^'\k - q- p) 6^'\p- q' - p') 



X 
X 



dri / dT2T{q,Ti,T2)'^^-^v{q)6^-^\q + q) 



Im + Pm) il'm' + P'm') ^hm' n' ^a' b' n' ^mni ^abn By {p') ■ (B.15) 

Using the three delta functions over the momenta the expression of Eq. flB.lSP becomes: 



HP{k,r) 



X 



(27r)3 

gaga' 
g' 



n 



d qv{q)By{k) / dri / (ir2 r(g, ri, 



km {km' qm') ^bm' n' ^a' b' n' ^mni ^abri' 



(B.16) 
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Eq. flB.16p can be also written, after some algebra, as 

H> '{k,T) = J d^q v{q) Bi[k) 

X r dn r dT2T{q,n,T2). (B.17) 
Jo Jo 

By appreciating that r(g, ri, = r(g, r2, Ti) the integrations over ti and T2 can be explicitly 
performed by recalling the elementary integral 

ry 



sin(6x) e ""^ dx 



7 ^ • 7 

cos by + - sm by 

b 



e-'^y. (B.18) 



The result is: 

{1 - e-^'^*h^[cos (gcsbr) + e(g, r) sin (gcsbr)]}^ 



JO Jo 



2q^cl[l + e'{q,T)]^ 

(B.19) 

where e{q,T) = q'^i^th/^^sh ^ 1- the limit €{q,T) <^ 1 and g^z/thT 7(q"r) — [Cgj^r^J/S and the 
higher order contributions 

k'^Biik) f d^q 



Hf(*.r) = (,..)„(,). 



2! 3 7 (27r)3 ' J (2 
1 A;2 5i(A^) f d^q /• d^p 

d^k 

^,l{k,T)v{k), (B.20) 
and so on can be resummed leading to an average magnetic field 

H,{K r) = e-^(^'-) B,{k), f{k, r) = ^^^jf^^An kl k^r' j . (B.21) 
The initial magnetic field is suppressed when f{k,r^) > 1, i.e. when 



J2(n + 1) ,,^/k \("+i)/2 
, >20 ^ , , (B.22) 

where dx is the (comoving) angular diameter distance to last scattering and knu denotes 
the diffusion scale in the non-Markovian approach proposed in this appendix. The results 
obtained here are consistent with a rough dimensional going, in short, as follows. Equation 
(IB.6P does not possess stationary solutions (see, e.g. [21]): the field is either amplified or 
dissipated. The typical diffusivity scale induced by large-scale curvature perturbations can 
be simply obtained by balancing the left and the right-hand sides of Eq. (IB.6P and by 
assuming that the typical amplitude of the velocity field is given by v^^ ~ V^n (Jot ^ I). 
In this case the typical diffusivity length is L ~ y/A-RT ~ y/A-jiH^^. 
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